CMB ANISOTROPIES: TOTAL ANGULAR MOMENTUM METHOD 
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A total angular momentum representation simplifies the radiation transport problem for tempera- 
ture and polarization anisotropy in the CMB. Scattering terms couple only the quadrupole moments 
of the distributions and each moment corresponds directly to the observable angular pattern on the 
sky. We develop and employ these techniques to study the general properties of anisotropy gen- 
eration from scalar, vector and tensor perturbations to the metric and the matter, both in the 
cosmological fluids and from any seed perturbations (e.g. defects) that may be present. The sim- 
pler, more transparent form and derivation of the Boltzmann equations brings out the geometric 
and model-independent aspects of temperature and polarization anisotropy formation. Large an- 
gle scalar polarization provides a robust means to distinguish between isocurvature and adiabatic 
models for structure formation in principle. Vector modes have the unique property that the CMB 
polarization is dominated by magnetic type parity at small angles (a factor of 6 in power compared 
with for the scalars and 8/13 for the tensors) and hence potentially distinguishable independent 
of the model for the seed. The tensor modes produce a different sign from the scalars and vectors 
for the temperature-polarization correlations at large angles. We explore conditions under which 
one perturbation type may dominate over the others including a detailed treatment of the photon- 
baryon fluid before recombination. 
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I. INTRODUCTION 



The cosmic microwave background (CMB) is fast becoming the premier laboratory for early universe and classical 
cosmology. With the flood of high quality data expected in the coming years, most notably from the new MAP Q| 
and Planck Surveyor [|| satellite missions, it is imperative that theoretical tools for their interpretation be developed. 
The corresponding techniques involved should be as physically transparent as possible so that the implications for 
cosmology will be readily apparent from the data. 

Toward this end, we reconsider the general problem of temperature and polarization anisotropy formation in the 
CMB. These anisotropics arise from gravitational perturbations which separate into scalar (compressional), vector 
(vortical), and tensor (gravity wave) modes. In previous treatments, the simple underlying geometrical distinctions and 
physical processes involved in their appearance as CMB anisotropics has been obscured by the choice of representation 
for the angular distribution of the CMB. In this paper, we systematically develop a new representation, the total 
angular momentum representation, which puts vector and tensor modes for the temperature and all polarization modes 
on an equal footing with the familiar scalar temperature modes. For polarization, this completes and substantially 
simplifies the ground-breaking work of PM|. Although we consider only flat geometries here for simplicity, the 
framework we establish allows for straightforward generalization to open geometries m-M unlike previous treatments. 

The central idea of this method is to employ only observable quantities, i.e. those which involve the total angular 
dependence of the temperature and polarization distributions. By applying this principle from beginning to end, we 
obtain a substantial simplification of the radiation transport problem underlying anisotropy formation. Scattering 
terms couple only the quadrupole moments of the temperature and polarization distributions. Each moment of the 
distribution corresponds to angular moments on the sky which allows a direct relation between the fundamental 
scattering and gravitational sources and the observable anisotropy through their integral solutions. 

We study the means by which gravitational perturbations of the scalar, vector, or tensor type, originating in either 
the cosmological fiuids or seed sources such as defects, form temperature and polarization anisotropies in the CMB. 
As is well established [|[|| , scalar perturbations generate only the so-called electric parity mode of the polarization. 
Here we show that conversely the ratio of magnetic to electric parity power is a factor of 6 for vectors, compared 
with 8/13 for tensors, independent of their source. Furthermore, the large angle limits of polarization must obey 
simple geometrical constraints for its amplitude that differ between scalars, vectors and tensors. The sense of the 
temperature-polarization cross correlation at large angle is also determined by geometric considerations which separate 
the scalars and vectors from the tensors [^- These constraints are important since large-angle polarization unlike 
large-angle temperature anisotropies allow one to see directly scales above the horizon at last scattering. Combined 
with causal constraints, they provide robust signatures of causal isocurvature models for structure formation such as 
cosmological defects. 

In §p we develop the formalism of the total angular momentum representation and lay the groundwork for the 
geometric interpretation of the radiation transport problem and its integral solutions. We further establish the 
relationship between scalars, vectors, and tensors and the orthogonal angular modes on the sphere. In §[II, we treat 



the radiation transport problem from first principles. The total angular momentum representation simplifies both the 
derivation and the form of the evolution equations for the radiation. We present the differential form of these equations, 
their integral solutions, and their geometric interpretation. In §|^ we specialize the treatment to the tight-coupling 
limit for the photon-baryon fluid before recombination and show how acoustic waves and vorticity are generated 
from metric perturbations and dissipated through the action of viscosity, polarization and heat conduction. In §^, 
we provide specific examples inspired by seeded models such as cosmological defects. We trace the full process that 



transfers seed fluctuations in the matter through metric perturbations to observable anisotropics in the temperature 
and polarization distributions. 

II. NORMAL MODES 

In this section, we introduce the total angular momentum representation for the normal modes of fluctuations in 
a flat universe that are used to describe the CMB temperature and polarization as well as the metric and matter 
fluctuations. This representation greatly simplifies the derivation and form of the evolution equations for fluctuations in 



III. In particular, the angular structure of modes corresponds directly to the angular distribution of the temperature 



and polarization, whereas the radial structure determines how distant sources contribute to this angular distribution. 
The new aspect of this approach is the isolation of the total angular dependence of the modes by combining 
the intrinsic angular structure with that of the plane-wave spatial dependence. This property implies that the 
normal modes correspond directly to angular structures on the sky as opposed to the commonly employed technique 
that isolates portions of the intrinsic angular dependence and hence a linear combination of observable modes [lOl . 
Elements of this approach can be found in earlier works (e.g. P,|7|,pl| for the temperature and for the scalar 
and tensor polarization) . We provide here a systematic study of this technique which also provides for a substantial 



simplification of the evolution equations and their integral solution in § [II C , including the terms involving the radiation 
transport of the CMB. We discuss in detail how the monopole, dipole and quadrupole sources that enter into the 
radiation transport problem project as anisotropics on the sky today. 

Readers not interested in the formal details may skip this section on first reading and simply note that the tem- 
perature and polarization distribution is decomposed into the modes Y^™ exp(ifc • x) and -1-2^'" exp(ifc • x) with 
m = 0,±1,±2 for scalar, vector and tensor metric perturbations respectively. In this representation, the geometric 
distinction between scalar, vector and tensor contributions to the anisotropics is clear as is the reason why they do not 
mix. Here the ±2^£™ ^^^ the spin-2 spherical harmonics ||l^ and were introduced to the study of CMB polarization 

by Q. The radial decompositions of the modes Y^j^,™' (kr) and ±2^™[4'" V^'") i ^/^l™ (^^)] (for £ = 2) isolate the 
total angular dependence by combining the intrinsic and plane wave angular momenta. 



A. Angular Modes 



In this section, we derive the basic properties of the angular modes of the temperature and polarization distributions 
that will be useful in §[11 to describe their evolution. In particular, the Clebsch-Gordan relation for the addition of 



angular momentum plays a central role in exposing the simplicity of the total angular momentum representation. 

A scalar, or spin-0 field on the sky such as the temperature can be decomposed into spherical harmonics Y^™. 
Likewise a spin-s field on the sky can be decomposed into the spin- weighted spherical harmonics ^Yf^ and a tensor 
constructed out of the basis vectors eg ± ie^, e^ |131. The basis for a spin-2 field such as the polarization is _|_2^/"^± 
|,| where 



M± 



(ee T ie^) (g) {eg =F ie^) 



(1) 



since it transforms under rotations as a 2 x 2 symmetric traceless tensor. This property is more easily seen through 



the relation to the Pauli matrices, M± = a^ T*fij in spherical coordinates 
in terms of rotation matricesn as |l2] 



The spin-s harmonics are expressed 
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TAB. 1: Quadrupole {£ = 2) harmonics for spin-0 and 2. 



The rotation matrix 2?!^ ,„((/), 0, V") = ^/I^T/i2^^^ ^Y^" {0 , (j))e~'''''^ represents rotations by the Euler angles (0, 0, ^). 
Since the spin-2 harmonics will be useful in the following sections, we give their explicit form in Table 1 for £ = 2; the 
higher i harmonics are related to the ordinary spherical harmonics as 
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By virtue of their relation to the rotation matrices, the spin harmonics satisfy: the compatibility relation with 
spherical harmonics, qY™ — Y™\ the conjugation relation ^Y^™* == (—1)™+* _gY^^™; the orthonormality relation, 
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the completeness relation. 
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the parity relation. 






the generalized addition relation, 

Y,[,^Yr {0', <!>')] [s^ypio,. 
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(6) 
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which follows from the group multiplication property of rotation matrices which relates a rotation from 



through 



the origin to {0, (p) with a direct rotation in terms of the Eulcr angles (a, /3, 7) defined in Fig. |l|; and the Clebsch- 
Gordan relation, 
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It is worthwhile to examine the implications of these properties. Note that the orthogonality and completeness 
relations Eqns. (||) and (||) do not extend to different spin states. Orthogonality between s = ±2 states is established 
by the Pauh basis of Eqn. (||) M^M_|_ = 1 and M^M-p = 0. The parity equation (||) tells us that the spin flips 
under a parity transformation so that unlike the s = spherical harmonics, the higher spin harmonics are not parity 
eigenstates. Orthonormal parity states can be constructed as 




FIG. 1. Addition theorem and scattering geometry. The addition theorem for spin-s harmonics Eqn. (|7|) is established by 
their relation to rotations Eqn. (EI) and by noting that a rotation from {6',4>') through the origin (pole) to {6,(f>) is equivalent 
to a direct rotation by the Euler angles (a,/?, 7). For the scattering problem of Eqn. (ES), these angles represent the rotation 
by a from the k — e^ frame to the scattering frame, by the scattering angle /3, and by 7 back into the k frame. 
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[,YrM+± _,YrM^], 



(9) 



which have "electric" (—1)^ and "magnetic" (—1)^+^ type parity for the (±) states respectively. We shall see in §111 C 
that the polarization evolution naturally separates into parit y eige nstates. The addition property will be useful in 

l)is 



relating the scattering angle to coordinates on the sphere in §|IIB[ Finally the Clebsch-Gordan relation Eqn 
central to the following discussion and will be used to deri ve the total angular monicntum representation in i 
evolution equations for angular moments of the radiation §[IIC. 
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B. Radial Modes 



We now complete the formalism needed to describe the temperature and polarization fields by adding a spatial 
dependence to the modes. By further separating the radial dependence of the modes, we gain insight on their full 
angular structure. This decomposition will be useful in constructing the formal integral solutions of the perturbation 
equations in § [II C . We begin with its derivation and then proceed to its geometric interpretation. 



1. Derivation 



The temperature and polarization distribution of the radiation is in general a function of both spatial position x 
and angle n defining the propagation direction. In flat space, we know that plane waves form a complete basis for the 
spatial dependence. Thus a spin-0 field like the temperature may be expanded in 




FIG. 2. Projection effects. A plane wave exp(iA;-a;) can be decomposed into jt{kr)Yg and hence carries an "orbital" angular 
dependence. A plane wave source at distance r thus contributes angularpower to I ~ kr aX 9 = 7r/2 but also to larger angles 
^ <C fcr at ^ = which is encapsulated into the structure of ji (see Fig. H). If the source has an intrinsic angular dependence, 
the distribution of power is altered. For an aligned dipole Y^ oc cos 9 ('figure 8's) power at 6 = 7r/2 or £ ^ kr is suppressed. 
These arguments are generalized for other intrinsic angular dependences in the text. 



GT = HY 
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17"(n)exp(ifc • x) , 



(10) 



where the normahzation is chosen to agree vifith the standard Legendre polynomial conventions for m ^ 0. Likewise 
a spin-2 field like the polarization may be expanded in 
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The plane wave itself also carries an angular dependence of course, 

exp(ifc • x) = ^(-i) V47r(2^ + l)jt{kr)Y^{n) , 



(11) 



(12) 



where 63 = fc and x = —rfi (see Fig. 0). The sign convention for the direction is opposite to direction on the sky to 
be in accord with the direction of propagation of the radiation to the observer. Thus the extra factor of (—1)^ comes 
from the parity relation Eqn. m). 

The separation of the mode functions into an intrinsic angular dependence and plane-wave spatial dependence is 
essentially a division into spin (gY^) and orbital (Y^) angular momentum. Since only the total angular dependence 
is observable, it is instructive to employ the Clebsch-Gordan relation of Eqn. (H) to add the angular momenta. In 
general this couples the states between \£ — i'\ and £ + £'. Correspondingly a state of definite total £ will correspond 
to a weighted sum of j\e-e'\ to jt+i' in its radial dependence. This can be reexpressed in terms of the je using the 
recursion relations of spherical Bessel functions. 



We can then rewrite 
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FIG. 3. Radial spin-0 (temperature) modes. The angular power in a plane wave (left panel, top) is modified due to the 
intrinsic angular structure of the source as discussed in the text. The left panel corresponds to the power in scalar {m = 0) 
monopole Gq, dipole G?, and quadrupole G% sources (top to bottom); the right panel to that in vector {m — 1) dipole G-^ and 
quadrupole Gj sources and a tensor (m — 2) quadrupole Gj source (top to bottom). Note the differences in how sharply 
peaked the power is at £ « kr and how fast power falls as € <S kr. The argument of the radial functions kr = 100 here. 



where the lowest {£',m) radial functions are 
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with primes representing derivatives with respect to the argument of the radial function x — kr. These modes are 
shown in Fig. 0. 

Similarly for the spin ±2 functions with ra > Q (see Fig. H), 



±2G^ = 5](-*)V4^(2^+l)[4™^(fc^) ± ^P^^\kr)] ^^Y^in) , 
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where 
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which corresponds to the £' = i,£±2 coupling and 
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FIG. 4. Radial spin-2 (polarization) modes. Displayed is the angular power in a plane-wave spin-2 source. The top panel 
shows that vector {m = 1, upper panel) sources are dominated by _B-parity contributions, whereas tensor (m — 2, lower panel) 
sources have comparable but less power in the B-parity. Note that the power is strongly peaked sX I = kr for the B-parity 
vectors and i5-parity tensors. The argument of the radial functions kr = 100 here. 
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which corresponds to the £' = ^ ± 1 coupling. The corresponding relation for negative m involves a reversal in sign of 
the /3-functions 



These functions are plotted in Fig. ^. Note that e\ ^ i\ is displayed in Fig 
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2. Interpretation 



The structure of these functions is readily apparent from geometrical considerations. A single plane wave contributes 
to a range of angular scales from (, k, kr B,t 9 = 7r/2 to larger angles i ^ kr as, ^r (0,7r), where k ■ h = cos9 (see 
Fig. mj. The power in ^ of a single plane wave shown in Fig. S(a) (top panel) drops to zero i > kr, has a concentration 
of power around £ = kr and an extended low amplitude tail to £ ^ kr. 

Now if the plane wave is multiplied by an intrinsic angular dependence, the projected power changes. The key to 
understanding this effect is to note that the intrinsic angular behavior is related to power in £ as 



e-^TT/2 ^=> (.f=ikr. (20) 

Thus factors of sin^ in the intrinsic angular dependence suppress power at £ <C fcr ("aHasing suppression"), whereas 
factors of cos0 suppress power at £ ~ fcr ("projection suppression"). Let us consider first a m = dipole contribution 



Y]° ex cos9 (see Fig 
£ fn kr (compare Fig. 3 



2). The cos 9 dependence suppresses power in j) at the peak in the plane- wave spectrum 
(a) top and middle panels). The remaining power is broadly distributed for i ^ kr. The same 
reasoning applies for Y2 quadrupole sources which have an intrinsic angular dependence of 3cos^0 — 1. Now the 
minimum falls at = cos~'^(l/V3) causing the double peaked form of the power in jj shown in Fig. |2(a) (bottom 
panel). This series can be continued to higher G^ and such techniques have been used in the free streaming limit for 
temperature anisotropics |11[ . 

Similarly, the structures of j^ i ii ^-nd i\ are apparent from the intrinsic angular dependences of the G\ , G\ 
and G\ sources, 

Y^ oc sin 6*6*"* , Y^ ex sin 9 cos 6*6"^ , ^3^ ex sin^ ^e^*"* , (21) 

respectively. The sin0 factors imply that as m increases, low (, power in the source decreases (compare Fig. @(a,b) 
top panels). G\ suffers a further suppression aX 9 = 7r/2 (£ « kr) from its cos0 factor. 

There are two interesting consequences of this behavior. The sharpness of the radial function around £ ~ kr 
quantifies how faithfully features in the /c-space spectrum are preserved in ^-space. If all else is equal, this faithfulness 
increases with |?7i| for GV^, due to aliasing suppression from sin™ 9. On the other hand, features in Gl^ji , j^ are washed 
out in comparison due projection suppression from the cos9 factor. 

Secondly, even if there are no contributions from long wavelength sources with k <C £/r, there will still be large 
angle anisotropics at £ <C fcr which scale as 

[£jf™']2(x^2+2|™|^ (22) 

This scaling puts an upper bound on how steeply the power can rise with £ that increases with \m\ and hence a lower 
bound on the amount of large relative to small angle power that decreases with \m\. 

The same arguments apply to the spin-2 functions with the added complication of the appearance of two radial 
functions eg and /?£. The addition of spin-2 angular momenta introduces a /3-contribution from e*""^ except for m — 0. 
For m = ±1, the /3-contribution strongly dominates over the e-contributions; whereas for m = ±2, e-contributions are 
slightly larger than /3-contributions (see Fig. ^. The ratios reach the asymptotic values of 

for fixed kr ^ 1. These considerations are closely related to the parity of the multipole expansion. Although the 
orbital angular momentum does not mix states of different spin, it does mix states of different parity since the plane 
wave itself does not have definite parity. A state with "electric" parity in the intrinsic angular dependence (see Eqn. 0) 
becomes 

2G^M+ + _2G?M_ = ^(-*)V47r(2^+ l){4"^[2>^rM+ + _2>^rM-] 

+tp^r\2YrM+~ _,YrM^]}. (24) 

Thus the addition of angular momentum of the plane wave generates "magnetic" _B-type parity of amplitude (3i out of 
an intrinsically "electric" i?-type source as well as _E-type parity of amplitude eg. Thus the behavior of the two radial 



functions has significant consequences for the polarization calculation in § [II C and implies that B-parity polarization 
is absent for scalars, dominant for vectors, and comparable to but slightly smaller than the E-parity for tensors. 

Now let us consider the low £ ^ kr tail of the spin-2 radial functions. Unlike the spin-0 projection, the spin-2 
projection allows increasingly more power at — > and/or tt, i.e. £ <^ kr, as \m\ increases (see Table 1 and note the 
factors of sin0). In this limit, the power in a plane wave fluctuation goes as 

[^e^'"^]^ ex f6-2|m|^ [^^M]2 ^ ^6-2|™|^ (25) 




FIG. 5. Spin-0 x Spin-2 (temperature x polarization) modes. Displayed is the cross angular power in plane wave spin-0 
and spin-2 sources. The top panel shows that a scalar monopole (m — 0) source correlates with a scalar spin-2 (polarization 
quadrupole) source whereas the tensor quadrupole (m — 2) anticorrelates with a tensor spin-2 source. Vector dipole (m — 1) 
sources oscillate in their correlation with vector spin-2 sources and contribute negligible once modes are superimposed. One 
must go to vector quadrupole sources (lower panel) for a strong correlation. The argument of the radial functions kr = 100 
here. 



Comparing these expressions with Eqn. ( |22| ) , we note that the spin-0 and spin-2 functions have an opposite dependence 
on m. The consequence is that the relative power in large vs. small angle polarization tends to decrease from the 
m = 2 tensors to the m = scalars. 

Finally it is interesting to consider the cross power between spin-0 and spin-2 sources because it will be used to 
represent the temperature-polarization cross correlation. Again interesting geometric effects can be identified (see 



,•(00) (0) 



correlates (Fig. ^, top panel solid line, positive definite); for to = 1, 



.(22 



Ji 



(ll)Jl) 



Fig. 0). For TO, = 0, the power in 

oscillates (short dashed line) and for m — 2, jf'^'ef anticorrelates (long dashed line, negative definite). The cross 

power involves only el' jj, due to the opposite pa rity o f the pf^ modes. 



These properties will become important in § [11 and [VE and translates into cross power contributions with opposite 
sign between the scalar monopole temperature cross polarization sources and tensor quadrupole temperature cross 
polarization sources |^. Vector dipole temperature and polarization sources do not contribute strongly to the cross 
power since correlations and anticorrelations in j^ e^ will cancel when modes are superimposed. The same is true 
of the scalar dipole temperature cross polarization j^ e\ as is apparent from Figs. and 0. The vector cross power 
is dominated by quadrupole temperature and polarization sources j\ ej, (Fig. ^ lower panel). 



C. Perturbation Classification 

As is well known (see e.g. HJlS]), a general symmetric tensor such as the metric and stress-energy perturbations 
can be separated into scalar, vector and tensor pieces through their coordinate transformation properties. We now 
review the properties of their normal modes so that they may be related to those of the radiation. We find that 
the TO = 0,±1,±2 modes of the radiation couple to the scalar, vector and tensor modes of the metric. Although 



10 



we consider flat geometries here, we preserve a covariant notation that ensures straightforward generaHzation to 
open geometries through the replacement of Sij with the curved three metric and ordinary derivatives with covariant 
derivatives HQ- 

1. Scalar Perturbations 

Scalar perturbations in a flat universe are represented by plane waves Q^'~'' — exp{ik ■ x) , which are the eigenfunctions 
of the Laplacian operator 

V2g(o)=-fc2g(o)^ (26) 

and their spatial derivatives. For example, vector and symmetric tensor quantities such as velocities and stresses 
based on scalar perturbations can be constructed as 

gf ) = -fc-iV,Q(°), Qf = [fc-^V,V, + i<5.,]Q(°). (27) 

Since V x Q^"-* = 0, velocity flelds based on scalar perturbations are irrotational. Notice that Q*^"^ — Gq, n^Q\ ~ G\ 
and rC"n?Q\; oc G\, where the coordinate system is defined by 63 = fc. From the orthogonality of the spherical 
harmonics, it follows that scalars generate only ttt, = fluctuations in the radiation. 

2,. Vector Perturbations 

Vector perturbations can be decomposed into harmonic modes Q^ of the Laplacian in the same manner as the 
scalars, 

V^gi^^) = -fc2Qp\ (28) 



which satisfy a divergenceless condition 



V'gf'^=0. (29) 



A velocity field based on vector perturbations thus represents vorticity, whereas scalar objects such as density per- 
turbations are entirely absent. The corresponding symmetric tensor is constructed out of derivatives as 



gr = -^(v.gr) + v,gr)). (3o) 



A convenient representation is 



Ql — p(ei ±ie2)iexp(z/c • x) . (31) 

v2 

Notice that n^Q\ = Gi^ and n^n^Qlj ex Gf^. Thus vector perturbations stimulate the m = ±1 modes in the 
radiation. 

3. Tensor Perturbations 
Tensor perturbations are represented by Laplacian eigenfunctions 

V2g(f)=-fc2g(±2)^ (32) 

which satisfy a transverse-traceless condition 

(5'^glf^ = v'glf' = 0, (33) 
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that forbids the construction of scalar and vector objects such as density and velocity fields. The modes take on an 
explicit representation of 

/ q 

Qii = ^\/ 0(^1 ±«e2)j (8) (ei ± ie2)j exp(ifc • x) . (34) 

■^ V o 

Notice that tl'tl^Qi- = G2 and thus tensors stimulate the m = ±2 modes in the radiation. 

In the following sections, we often only explicitly show the positive m value with the understanding that its opposite 
takes on the same form except where otherwise noted (i.e. in the _B-typc polarization where a sign reversal occurs). 

III. PERTURBATION EVOLUTION 

We discuss here the evolution of perturbations in the normal modes of §0. We first review the decomposition of 
perturbations in the metric and stress-energy tensor into scalar, vector and tensor types (§|[II A). We further divide 



the stress-energy tensor into fluid contributions, applicable to the usual particle species, and seed perturbations, 
applicable to cosmological defect models. We then employ the techniques developed in §0 to obtain a new, simpler 
derivation and form of the radiation transport of the CMB under Thomson scattering, including polarization (§111 B|), 



than that obtained first by [h6|. The complete evolution equations, given in §[IIC, are again substantially simpler 
in form than those of prior works where they overlap 0fl,nn,ni| and treats the case of vector perturbations. Finally 



m 



[II D, we derive the formal integral solutions through the use of the radial functions of §|IH and discuss their 



geometric interpretation. These solutions encapsulate many of the important results. 

A. Perturbations 

1. Metric Tensor 
The ultimate source of CMB anisotropics is the gravitational redshift induced by the metric fluctuation h^^, 

gi^ii^ = a^iViJ,!^ + hf_,„) , (35) 

where the zeroth component represents conformal time drj — dt/a and, in the flat universe considered here, 77^^ is the 
Minkowski metric. The metric perturbation can be further broken up into the normal modes of scalar, vector and 



tensor types as in ^I C . Scalar and vector modes exhibit gauge freedom which is fixed by an explicit choice of the 
coordinates that relate the perturbation to the background. For the scalars, we choose the Newtonian gauge (see e.g. 

/loo = 2*Q('", % = 2$Q(")% , (36) 

where the metric is shear free. For the vectors, we choose 



and the tensors 



ho^^-VQ^^\ (37) 



/ly = 2HQf^. (38) 



Note that tensor fluctuations do not exhibit gauge freedom of this type. 

2. Stress Energy Tensor 

The stress energy tensor can be broken up into fluid (/) contributions and seed (s) contributions (see e.g. |lq| ). The 
latter is distinguished by the fact that the net effect can be viewed as a perturbation to the background. Specifically 
Tfj,,u = T^y -\- ST^i, where T^g = —pf, 7"° = Tg' = and T*^ = Pf^^j is given by the fluid alone. The fluctuations can 



be decomposed into the normal modes of §IIC as 
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for the scalar components, 



for the vector components, and 



for the tensor components. 



5T", = -[pfSf+p,]Q("\ 

ST\ = -[{Pf+Pf)vf^+Ps]Q'^'^\ 
5T\ = [{Pf+Pf){vf^-V)+v^^]Qf\ 



5r,. = [p,4^) + 7rf]g(^K 



(39) 



(40) 



(41) 



B. Radiation Transport 

1. Stokes Parameters 

The Boltzmann equation for the CMB describes the transport of the photons under Thomson scattering by the 
electrons. The radiation is described by the intensity matrix: the time average of the electric field tensor E*Ej over 
a time long compared to the frequency of the light or equivalently as the components of the photon density matrix 
(see |19|] for reviews). For radiation propagating radially -E _L e^, so that the intensity matrix exists on the eg (E) e^ 
subspace. The matrix can further be decomposed in terms of the 2x2 Pauli matrices ai and the unit matrix 1 on 
this subspace. 

For our purposes, it is convenient to describe the polarization in temperature fluctuation units rather than intensity, 
where the analogous matrix becomes, 



T^ei + Q(T3 + U(Ji + V(J2- 



(42) 



<d — Tr(Tl)/2 — AT/T is the temperature perturbation summed over polarization states. Since Q — Tr(Tcr3)/2, 
it is the difference in temperature fluctuations polarized in the eg and e^ directions. Similarly U = Tr(TCTi)/2 is 
the difference along axes rotated by 45°, (eg ± e0)/\/2, and V = Tr(Tcr2)/2 that between {eg ± ie^)/^/2. Q and 
U thus represent linearly polarized light in the north/south-east/west and northeast/southwest-northwest/southeast 
directions on the sphere respectively. V represents circularly polarized light (in this section only, not to be confused 
with vector metric perturbations). 

Under a counterclockwise rotation of the axes through an angle ip the intensity T transforms as T' — RTR~^. O 
and V remain distinct while Q and U transform into one another. Since the Pauli matrices transform as 173 ± ia'i = 



^^2iip 



((73 ± JCTi) a more convenient description is 



61 + y<T2 + {Q + ^C/)M+ + (Q - ^C/)M_ 



(43) 



where recall that M± — (0-3 =F iai)/2 (see Eqn. ||), so that Q±iU transforms into itself under rotation. Thus Eqn. (g) 
implies that Q ± ill should be decomposed into s = ±2 spin harmonics [^|4|. 

Since circular polarization cannot be generated by Thomson scattering alone, we shall hereafter ignore V. It is then 
convenient to reexpress the matrix as a vector 



T=ie,Q + iU,Q-iU). 



(44) 



The Boltzmann equation describes the evolution of the vector T under the Thomson coUisional term C[T] and 
gravitational redshifts in a perturbed metric G[hfj_^] 
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7 O 

— f (t], f , n) = —f + n'V^f = C[f] + G[h^.] , (45) 

dr] or] 

where we have used the fact that Xi = rii and that in a flat universe photons propagate in straight hnes h — Q. We 
shaU now evaluate the Thomson scattering and gravitational redshift terms. 

2. Scattering Matrix 

The calculation of Thomson scattering including polarization was first performed by Chandrasekhar [|l6| ; here we 
show a much simpler derivation employing the spin harmonics. The Thomson differential scattering cross section 
depends on angle as |e' • ep where e' and e are the incoming and outgoing polarization vectors respectively in the 
electron rest frame. Radiation polarized perpendicular to the scattering plane scatters isotropically, while that in the 
scattering plane picks up a factor of cos^ /3 where fi is the scattering angle. If the radiation has different intensities or 
temperatures at right angles, the radiation scattered into a given angle will be linearly polarized. 

Now let us evaluate the scattering term explicitly. The angular dependence of the scattering gives 




cos2/3 


\ 


/e,! 


1 





e^ 





cos/3 / 


\ u 



(46) 



where the U transformation follows from its definition in terms of the difference in intensities polarized ±45° from 



the scattering plane. With the relations 6 = 8|| 
T representation of Eqn. (44) becomes ,g 



8^ and Q ± it/ = 6|| — 0^ ± iU , the angular dependence in the 



f 



4 



cos^ (3 + 1 


-\ sm^ fi 


-isin^/? 


-isin^/? 


i(cos/3 + l)2 


i(cos/3-l)2 


-\siT? (3 


i (cos /3- 1)2 


\{cos(i+lf ^ 



T, 



(47) 



where the overall normalization is fixed by photon conservation in the scattering. To relate these scattering frame 
quantities to those in the frame defined by fc = 63, we must first perform a rotation from the k frame to the 
scattering frame. The geometry is displayed in Fig. H^, where the angle a separates the scattering plane from the 
meridian plane at {6',4>') spanned by e^ and eg. After scattering, we rotate by the angle between the scattering 
plane and the meridian plane at (6*, (j)) to return to the k frame. Eqn. ( p3| ) tells us these rotations transform T as 
R(V')f = diag(l, e^*'^, e-^**^)!^ . The net result is thus expressed as 



R(7)S(/3)R(-a) = \^j^ 



^ Y§{p,a) + 2VlYS{P,a) 


~^Y~^{P,a) 


-^Yi{(i,a) ^ 


-V&2Y§{P,a)e-^''^ 


Z^Y^-^{l3,a)e-^-^ 


3 2i^2^(/9,a)e-2n 


^ -V^_^Y°{(i,a)e^''^ 


i_,Y^-\t3,a)e'^^ 


2._^Yi{l3,a)e^''^ ^ 



(48) 



where we have employed the explict spin-2, £ = 2 forms in Tab. 1. Integrating over incoming angles, we obtain the 
collision term in the electron rest frame 



dD.' 

C[T]rcst = -Tf{n) + f / — R(7)S(/3)R(-a)T(f^') , 



(49) 



where the two terms on the rhs account for scattering out of and into a given angle respectively. Here the differential 
optical depth r — n^cTTa sets the collision rate in conformal time with Uf, as the free electron density and ctt as the 
Thomson cross section. 



Chandrasekhar employs a different sign convention for U 



-U. 
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The transformation from the electron rest frame into the background frame yields a Doppler shift n ■ vb in the 
temperature of the scattered radiation. With the help of the generalized addition relation for the harmonics Eqn. (|^), 
the full collision term can be written as 



C[f] 



--^^/-i/--'->'' 



(50) 



The vector / describes the isotropization of distribution in the electron rest frame and is given by 

dn' 



i{n) = T{n) 



An 



e' + n-VB, 0,0 



(51) 



The matrix p('") encapsulates the anisotropic nature of Thomson scattering and shows that as expected polarization 
is generated through quadrupole anisotropics in the temperature and vice versa 



p(™) 



/ 



V 






3 \r'ml \rm 
22^2 J^2 



•2-'2 2^2 



'-'2^2 -2-'^2 



3 \^mf \^m 
2-2^2 ^2 



\ 



'-2-'2 2-^2 



'-'-2^2 -2^2 



(52) 



/ 



where Y™' = Y^™*(ri') and ^i^™' = s^'"*(^') ^^d the unprimed harmonics are with respect t o fl. Th ese m 
0, ±1, ±2 components correspond to the scalar, vector and tensor scattering terms as discussed in §11 C and IIIC. 



3. Gravitational Redshift 

In a perturbed metric, gravitational interactions alter the temperature perturbation O. The redshift properties may 
be formally derived by employing the equation of motion for the photon energy p = —u^^p^ where u'^ is the 4-velocity 
of an observer at rest in the background frame and p'^ is the photon 4-momcntum. The Euler-Lagrange equations of 
motion for the photon and the requirement that \u'^\ = 1 result in 



P 



1 



1 



- = -n'n-'hij - n^hoi - -n^'Vihoo , 



(53) 



which differs from 1 20 ^ since we take h to be the photon propagation direction rather than the viewing direction of 
the observer. The first term is the cosmological reshift due to the expansion of the spatial metric; it does not affect 
temperature perturbations ST/T. The second term has a similar origin and is due to stretching of the spatial metric. 
The third and fourth term are the frame dragging and time dilation effects. 
Since gravitational redshift affects the different polarization states alike, 



G[V] = 



^kj + n'hoi + -n'W^hoo ,0 ,0 



—n n 
2 



(54) 



in the T basis. We now explicitly evaluate the Boltzmann equation for scalar, vector, and tensor metric fluctuations 
of Eqns. (|6l)-(||). 



C. Evolution Equations 



In this section, we derive the complete set of evolution equations for the temperature and polarization distribution 
in the scalar, vector, tensor decomposition of metric fluctuations. Though the scalar and tensor fluid results can be 
found elsewhere in the literature in a different form (see e.g. |11 1Q|), the total angular momentum representation 
substantially simplifies the form and aids in the interpretation of the results. The vector derivation is new to this 
work. 
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1. Angular Moments and Power 
The temperature and polarization fluctuations are expanded into the normal modes defined in ^I B J 

•J \ J D -wi — _9 



^ > e m=-2 



(55) 
{Q±^U){v,x,n) = j ^Y. E (i^l™^ ± ^<^) ±20^- 



A comparison with Eqn. (g) and (|4^ ) shows that E^^ and B^ represent polarization with electric (—1)^ and 
magnetic (— l)^"*"^ type parities respectively ^^- Because the temperature B^ has electric type parity, only E^'' 

couples directly to the temperature in the scattering sources. Note that BjT' and _E)™ represent polarizations with 
Q and U interchanged and thus represent polarization patterns rotated by 45°. A simple example is given by the 
m = modes. In the fc-frame, E^ represents a pure Q, or north/south-east/ west, polarization field whose amplitude 
depends on 0, e.g. sin^ 9 for (i ~ 2. B^ represents a pure t/, or northwest/southeast -northeast/southwest, polarization 
with the same dependence. 

The power spectra of temperature and polarization anisotropies today are defined as, e.g. Cf^ = (|afmP) for 
= ^ agrriYf™' with the average being over the (2£ + 1) 77i-values. Recalling the normalization of the mode functions 



from Eqn. (nfl) and (11), we obtain 

2 



•^ m=-2 

where X takes on the values 8, E and B. There is no cross correlation Cf^ or Cf^ due to parity [see Eqns. (o) 
and (|4|)]. We also employ the notation C^ for the m contributions individually. Note that B^ ' — here due 

to azimuthal symmetry in the transport problem so that Cg =0. 

As we shall now show, the m ~ 0,±1,±2 modes are stimulated by scalar, vector and tensor perturbations in the 
metric. The orthogonality of the spherical harmonics assures us that these modes are independent, and we now discuss 
the contributions separately. 

2. Free Streaming 

As the radiation free streams, gradients in the distribution produce anisotropies. For example, as photons from 
different temperature regions intersect on their trajectories, the temperature difference is reflected in the angular 
distribution. This effect is represented in the Boltzmann equation (|45| ) gradient term, 

_, _^ Air 

h-V ~^ih-k^i\ —kY^ . (57) 

V ^ 

which multiplies the intrinsic angular dependence of the temperature and polarization distributions, Yp and ±2^™ 
respectively, from the expansion Eqn. (pq ) and the angular basis of Eqns. (M) and dl^). Free streaming thus involves 
the Clebsch-Gordan relation of Eqn. (|8|) 

V S'^i^^'^^ ^' ./iWTTW^^ ' iie+l)^'''' '^ ^i2£+l)i2£ + 3)^' '^'^ ^ ^ 

which couples the i" to £ ± 1 moments of the distribution. Here the coupling coefficient is 



^Our conventions differ from i] as {2£ + l)Al^f^ = 49^ '^'/(Stt)^'^^ and similarly for A^^^] with 9^'^' -^ -Ef'^\~B[ 



and so C^^ = —C^ ' but with other power spectra the same. 
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,K™ = y(l2-r7^^2y(£2Tr^2)/^ 



(59) 



As we shall now see, the result of this streaming effect is an infinite hierarchy of coupled ^-moments that passes power 
from sources at low multipoles up the £-chain as time progresses. 



3. Boltzmann Equations 



The explicit form of the Boltzmann equations for the temperature and polarization follows directly from the Clebsch- 
Gordan relation of Eqn. (58). For the temperature (s = 0), 



el™^ 



(2£- 1)^^-1 (2^ + 3)^+1 



re 



(m) 



S, 



(m) 



{£> m). 



(60) 



The term in the square brackets is the free streaming effect that couples the ^-modes and tells us that in the absence 
of scattering power is transferred down the hierarchy when k-q ^ 1. This transferral merely represents geometrical 
projection of fluctuations on the scale corresponding to k at distance 77 which subtends an angle given hy d ^ krj. The 
main effect of scattering comes through the f0^™ term and implies an exponential suppression of anisotropics with 
optical depth in the absence of sources. The source S*)" accounts for the gravitational and residual scattering effects, 



:(o) 



fe 



ft(0) 



$. 



.(0) 



TV 



S\ 



(0) 



fc*. 



:.(0) 



fP(0) 



4^) = rPW 



(61) 



H . 



The presence of 0q represents the fact that an isotropic temperature fluctuation is not destroyed by scattering. 
The Doppler effect enters the dipole [H. ~ \) equation through the baryon velocity Vg term. Finally the anisotropic 
nature of Compton scattering is expressed through 



1 
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p(-^) = ±. ei^^-Ves 



(m) 



(62) 



and involves the quadrupole moments of the temperature and i<^-polarization distribution only. 

The polarization evolution follows a similar pattern for ^ > 2, m > from Eqn. (|58| ) with s — ±2q 



S'f''^ = k 



M""^ = k 






2to 



(2^-1) 



^ r(™) 



[21 ~1) 



2m 

£{£ + !) 



B 



(m) 



-E, 



(m) 



_2JH+2_p(m) 
(2^ + 3) ^+1 






(2^ + 3) 



■-[E^P^ + \/6P(™)5,,2 



tB 



(m) 



(63) 
(64) 



Notice that the source of polarizati on p (™) enters only in the £'-mode quadrupole due to the opposite parity of 02 
and i?2. However, as discussed in §[IB, free streaming or projection couples the two parities except for the m = Q 



scalars. Thus S); = by geometry regardless of the source. It is unnecessary to solve separately for the m = — |m| 
relations since they satisfy the same equations and solutions with B 



(.-\m\) _ 



B}™'" and all other quantities equal. 



To complete these equations, we need to express the evolution of the metric sources ($, ^f", V, H). It is to this subject 
we now turn. 



*The expressions above were all derived assuming a flat spatial geometry. In this formalism, including the effects of spatial 
curvature is straightforward: the f ± 1 terms in the hierarchy are multiplied by factors of [1 — {£'^ — m— VjK/k^]^'^ 10101' where 
the curvature is K — ~Hq{1 — fitot). These factors account for geodesic deviation and alter the transfer of power through the 
hierarchy. A full treatment of such effects will be provided in |8| . 
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4- Scalar Einstein Equations 
The Einstein equations Gf^u = SvrGT^^ express the metric evolution in terms of the matter sources. With the form 



of the scalar metric and stress energy tensor given in Eqns. (36) and (|39[), the "Poisson" equations become 



fc2$ = 47rG'a2 
where the corresponding matter evolution is given by covariant conservation of the stress energy tensor T^^, 



(65) 



(0) 



5f ^ -(1 + wf){kvy' + 3$) - 3-6wf 



a 



drj 



„(o) 



{i + wf)vy' ^(i + wf) 



fc*- -(l-3w/)' 



,(0) 



+ Wfk{Spf/pf- -TTf), 



(66) 



for the fluid part, where Wf = pf/pf. These equations express energy and momentum density conservation respec- 
tively. They remain true for each fluid individually in the absence of momentum exchange. Note that for the photons 



_ 40(0) „(o) _ «(o) 



JO) _ 12o(0) 



, w-y =01 and TTj ' — -y02 • Massless neutrinos obey Eqn. ( |60| ) without the Thomson coupling term. 
Momentum exchange between the baryons and photons due to Thomson scattering follows by noting that for a 
given velocity perturbation the momentum density ratio between the two fluids is 

PB + PB 3pB 



P-r + Pj W't 
A comparison with photon Euler equation (pCJ ) (with t — 1, m — Q) gives the baryon equations as 



(67) 



-kv 



(0) 
B 



•,(0) 






3$, 



For a seed source, the conservation equations become 

Ps = -i-{Ps+Ps)-kvf^ 
a 

2 



^(o)^_4«.(o)+fc(p^_^,(o))^ 
a 6 



since the metric fluctuations produce higher order terms. 



(68) 



(69) 



5. Vector Einstein Equations 

The vector metric source evolution is similarly constructed from a "Poisson" equation 

V + 2-V = -8^Ga2(p/^(/) + 7r(i')/fc , 



and the momentum conservation equation for the stress-energy tensor or Euler equation 



v^P = V 



„2^«^,(l) 



^1^ ^/ -(1) 



(i-3c)-(.^^-y)-^fc-^.y 



^,W=-4^;(l)-ifc^(l), 



(70) 



(71) 



where recall ci — Pf/pf is the sound speed. Again, the first of these equations remains true for each fluid individually 

save for momentum exchange terms. For the photons v^ — Q\ and Triy — Iv302 ■ Thus with the photon Euler 
equation (pO[) (with £ = 1, rn = 1), the full baryon equation becomes 



•,(1) 



«^„,(i) 



L^fiW_„(i)^ 



a H 



(72) 



see Eqn. (pSh for details. 
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6. Tensor Einstein Equations 
The Einstein equations tell us that the tensor metric source is governed by 

H + 2-H + en = SnGa^p^Trf + ^(2)] _ (73) 



a 



where note that the photon contribution is tt^ — IO2 • 



D. Integral Solutions 

The Boltzmann equations have formal integral solutions that are simple to write down by considering the properties 



of source projection from § II B . The hierarchy equations for the temperature distribution Eqn. (Kd) merely express 



2^+1 


Jo ^ 




riv) s / f{r,')drj' 

Jrj 



the projection of the various plane wave temperature sources S^'G^ on the sky today (see Eqn. (|61|)). From the 
angular decomposition of G™ in Eqn. (|l4| ) , the integral solution immediately follows 

\Hvo~v))- (74) 

Here 

(75) 

is the optical depth between rj and the present. The combination fe~'^ is the visibility function and expresses the 
probability that a photon last scattered between drj of rj and hence is sharply peaked at the last scattering epoch. 
Similarly, the polarization solutions follow from the radial decomposition of the 

- V6fp(™) [2G2M+ + _2G^M_] (76) 

source. From Eqn. (p4), the solutions, 

^^^^^ = -V-^fJ d^re-p(™)(^)4")(fc(^o - V)) , 

^^Q^ = -Ve^" dr^re-^P^'-\r^)P^r\Km " v)) ■ (77) 

immediately follow as well. 

Thus the structures of y^ *" , e^™ , and /S)™ shown in Figs. || and Q directly reflects the angular power of the 

sources S^, and p(™). There are several general results that can be read off the radial functions. Regardless of the 
source behavior in fc, the _B-parity polarization for scalars vanishes, dominates by a factor of 6 over the electric parity 
at ^ ^ 2 for the vectors, and is reduced by a factor of 8/13 for the tensors at £ ^ 2 [see Eqn. (p3|)]. 
Furthermore, the power spectra in I can rise no faster than 

^2^ee(m) ^ ^2+2|m| g2fjEE{m) ^ ^6-2\m\ 

^^Cf ^(™) oc f 6-2|™l , ^^Cf ^("^ oc i\ (78) 

due to the aliasing of plane- wave power to £ <t^ k{rja — rj) (see Eqn. |25) which leads to interesting constraints on scalar 



temperature fluctuations [ |22[ and polarization fluctuations (see §VC). 

Features in fc-space in the £ — \m\ moment at fixed time are increasingly well preserved in ^-space as |?Ti| increases, 
but may be washed out if the source is not well localized in time. Only sources involving the visibility function re^"^ 
are required to be well localized at last scattering. However even features in such sources will be washed out if they 
occur in the £ = |to| + 1 moment, such as the scalar dipole and the vector quadrupole (see Fig. 0). Similarly features 
in the vector E and tensor B modes are washed out. 

The geometric properties of the temperature-polarization cross power spectrum Cf^ can also be read off the integral 
solutions. It is first instructive however to rewrite the integral solutions as {£ > 2) 
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9r(r/o,fc) 
2£+l 

2^+1 

Qf\vo,k) 
2^+1 



dije 



Vo 



di] e ^ 



Vo 



drje '^ 



(fe(°) + i^ + i>- $)jf °) + f.i°)j(^") + rP(")jf °) 



r(«i^' - V)^^^ + {rPi^^ + ^^kV)jr' 



f F(2) - H 



.•(22) 



(79) 



(0) 



^ 



where we have integrated the scalar and vector equations by parts noting that de '^ /drj — re '^ . Notice that 

acts as an effective temperature by accounting for the gravitational redshift from the potential wells at last scattering. 



We shall see in §IV that v^-^' « y at last scattering which suppresses the first term in the vector equation. Moreover, 



g-i 

,•(21) 



as discussed in §11 B and shown in Fig. H, the vector dipole terms (j) ) do not correlate well with the polarization 



(11) 



(e^ ) whereas the quadrupole terms iSt^') '^'^■ 

The cross power spectrum contains two pieces: the relation between the temperature and polarization sources of, 
and p('") respectively and the differences in their projection as anisotropics on the sky. The latter is independent 
of the model and provides interesting consequences in conjunction with tight coupling and causal constraints on the 
sources. In particular, the sign of the correlation is determined by pl| ] 

sgn[Cf^(°)] = -sgn[p(o)(e(") + v|/)], 
sgn[Cf ^(')] = -sgn[p(i)(\/3fp(i) + kV)\ , 
sgn[Cf ^*'^] = sgn[p(2) (fp(2) - ii) 



(80) 



where the sources are evaluated at last scattering and we have assumed that |8q + ^| 3> IP*-"-*! as is the case for 
standard recombination (see j|rv|). The scalar Doppler effect couples only weakly to the polarization due to differences 



in the projection (see §[IB). The important aspect is that relative to the sources, the tensor cross spectrum has an 
opposite sign due to the projection (see Fig. ^. 

These integral solutions are also useful in calculations. For example, they may be employed with approximate 
solutions to the sources in the tight coupling regime to gain physical insight on anisotropy formation (see §[V and 
|g2|,g3[). Seljak & Zaldarriaga |gj| have obtained exact solutions through numerically tracking the evolution of the 
source by solving the truncated Boltzmann hierarchy equations. Our expression agree with [plHJ24| where they overlap. 



IV. PHOTON-BARYON FLUID 



Before recombination, Thomson scattering between the photons and electrons and coulomb interactions between 
the electrons and baryons were sufhciently rapid that the photon-baryon system behaves as a single tightly coupled 
fluid. Formally, one expands the evolution equations in powers of the Thomson mean free path over the wavelength 
and horizon scale. Here we briefly review well known results for the scalars (see e.g. 12^J2^) to show how vector or 
vorticity perturbations differ in their behavior (§[VA). In particular, the lack of pressure support for the vorticity 
changes the relation between the CMB and metric fluctuations. We then study the h igher order effects of shear 
viscosity and polarization generation from scalar, vector and tensor perturbations (§ IVB ). We identify signatures in 
the temperature-polarization power spectra that can help separate the types of perturbations. Entropy generation 
a nd he at conduction only occur for the scalars (§ [VC ) and leads to differences in the dissipation rate for fluctuations 

(1x1). 



A. Compression and Vorticity 



For the {m = 0) scalars, the well-known result of expanding the Boltzmann equations 
baryon Euler equation (pSh is 



for £ = 0, 1 and the 
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which represent the photon fluid continuity and Euler equations and gives the baryon fluid quantities directly as 

to lowest order. Here mcs = I + R where recall that R is the baryon-photon momcntuni density ratio. We have 
dropped the viscosity term 02 = 0{k/T)Q[ (see § I V B ) . The effect of the baryons is to introduce a Compton drag 



term that slows the oscillation and enhances infall into gravitational potential wells ^. That these equations describe 
forced acoustic oscillations in the fluid is clear when we rewrite the equations as 

(meffe[,"^y+ ye^°) = -ymeff* - (meff4>y, (83) 

whose solution in the absence of metric fluctuations is 

6q — Ato"^' cos{ks + (f)) , 

e^°'> ^ V3AmJ^'^ sm{ks + 0) , (84) 

where s = J c-^sdri = J(3'm,cs)~^'^dri is the sound horizon, A is a constant amplitude and (f> is a. constant phase shift. 
In the presence of potential perturbations, the redshift a photon experiences climbing out of a potential well makes 
the effective temperature 6g + ^ (see Eqn. ^, which satisfies 

Kff(0^"^ + n+ y(0[,''^ + *) = -y^* + Kff(* - $)]•, (85) 

and shows that the effective force on the oscillator is due to baryon drag R'i? and differential gravitational redshifts 
from the time dependence of the metric. As seen in Eqn. (|7^) and (]8^), the effective temperature at last scattering 



forms the main contribution at last scattering with the Doppler effect v^' — Q]^ ' playing a secondary role for nics > 1. 
Furthermore, because of the nature of the monopole versus dipole projection, features in £ space are mainly created 
by the effective temperature (see Fig. || and § [II D| ) . 



If i? ^ 1, then one expects contributions of 0(^ — $)/fc^ to the oscillations in 8q + ^ in addition to the initial 
fluctuations. These acoustic contributions should be compared with the ©(A^ — A$) contributions from gravitational 
redshifts in a time dependent metric after last scattering. The stimulation of oscillations at fcry ^ 1 thus either requires 
large or rapidly varying metric fluctuations. In the case of the former, acoustic oscillations would be small compared 
to gravitational redshift contributions. 

Vector perturbations on the other hand lack pressure support and cannot generate acoustic or compressional waves. 
The tight coupling expansion of the photon (£ — l,m — 1) and baryon Euler equations (^ and (JT^ ) leads to 

Kff(0l'^ - "^)] ■ =0, (86) 

and Vg = ©1 ■ Thus the vorticity in the photon baryon fluid is of equal amplitude to the vector metric perturbation. 
In the absence of sources, it is constant in a photon-dominated fluid and decays as a^^ with the expansion in a baryon- 
dominated fluid. In the presence of sources, the solution is 

&['\v, k) = F(r7, k) + — [e(i)(0, k) - ViQ, k)] , (87) 

TOoff 

so that the photon dipole tracks the evolution of the metric fluctuation. With Vg = Q\' vn Eqn. (|6l|), vorticity leads 



to a Doppler effect in the CMB of magnitude on order the vector metric fluctuation at last scattering V in contrast 
to scalar acoustic effects which depend on the time rate of change of the metric. 
In turn the vector metric depends on the vector anisotropic stress of the matter as 



T/(?7,, k) = -87rGa;2 A drja'^ipfTrf'^ + ■K^P)lk. 



ex a^^ and decays with the expansion. They are thus generally negligible if the universe 
contains only the usual fluids. Only seeded models such as cosmological defects may have their contributions to the 
CMB anisotropy dominated by vector modes. However even though the vector to scalar fluid contribution to the 
anisotropy for seeded models is of order k'^V/{^ — $) and may be large, the vector to scalar gravitational redshift 
contributions, of order V/{'^ — $) is not necessarily large. Furthermore from the integral solution for the vectors 
Eqn. ( [7^ ) and the tight coupling approximation Eqn. (pTj), the fluid effects tend to cancel part of the gravitational 
effect. 
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FIG. 6. Power spectra for the standard cold dark matter model (scale invariant scalar adiabatic initial conditions with 
flo = 1, h = 0.5 and Q,Bh^ = 0.0125). Notice that B-polarization is absent, _E-polarization scales as £■* at large angles and the 
cross correlation (QE) is negative at large angles and reflects the acoustic oscillati ons a t small angles. In particular the phase 
of the EE and &E acoustic peaks is set by the temperature oscillations 00 (see §|[ V B ) . 



B. Viscosity and Polarization 



Anisotropic stress represents shear viscosity in the fluid and is generated as tight coupling breaks down on small 
scales where the photon diffusion length is comparable to the wavelength. For the photons, anisotropic stress is related 
to the quadrupole moments of the distribution 02 which is in turn coupled to the E-parity polarization £'2 • The 
zeroth order expansion of the polarization {£ — 2) equations (Eqn. B4) gives 
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^^2"'^ The quadrupole {i — 2) component of the temperature hierarchy (Eqn. 60) then becomes to lowest 
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for scalars and vectors. In the tight coupling limit, the scalar and vector sources of polarization traces the structure 
of the photon-baryon fluid velocity. For the tensors. 
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Combining Eqns. (p3) and (BCj), we see that polarization fluctuations are generally suppressed with respect to metric or 
temperature fluctuations. They are proportional to the quadrupole moments in the temperature which are suppressed 
by scattering. Only as the optical depth decreases can polarization be generated by scattering. Yet then the fraction 
of photons affected also decreases. In the standard cold dark matter model, the polarization is less than 5% of the 
temperature anisotropy at its peak (see Fig. M). 
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These scaling relations between the metric and anisotropic scattering sources of the temperature and polarization 
are important for understanding the large angle behavior of the polarization and temperature polarization cross 
spectrum. Here last scattering is effectively instantaneous compared with the scale of the perturbation and the tight 
coupling remains a good approximation through last scattering. 

For the scalars, the Euler equation ( pi| ) may be used to express the scalar velocity and hence the polarization in 
terms of the effective temperature, 
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mcs^)dr] . 
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Since rric 



1, Oj^ has the same sign as 0q 



conditions. It then follows that p(°) is also of the same sign and is of order 



^P before 0q + ^ itself can change signs, assuming reasonable initial 



(93) 



which is strongly suppressed for krj <^1. The definite sign leads to a definite prediction for the sign of the temperature 
polarization cross correlation on large angles. 
For the vectors 
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(94) 



and is both suppressed and has a definite sign in relation to the metric fluctuation. The tensor relation to the metric 
is given in Eqn. (K^). In fact, in all three cases the dominant source of temperature perturbations has the same sign 
as the anisotropic scattering source p(™). From Eqn. (80), differences in the sense of the cross correlation between 
temperature and polarization thus arise only due to geometric reasons in the projection of the sources (see Fig. 0). 
On angles larger than the horizon at last scattering, the scalar and vector Cf^ is negative whereas the tensor cross 
power is positive |P,pl|. 

On smaller scales, the scalar polarization follows the velocity in the tight coupling regime. It is instructive to recall 
the solutions for the acoustic oscillations from Eqn. (|8J). The velocity oscillates it/ 2 out of phase with the temperature 
and hence the £'-polarization acoustic peaks will be out of phase with the temperature peaks (see Fig. g). The cross 
correlation oscillates as cos(fcs + 0) sin(fcs + 0) and hence has twice the frequency. Thus between peaks of the 
polarization and temperature power spectra (which represents both peaks and troughs of the temperature amplitude) 
the cross correlation peaks. The structure of the cross correlation can be used to measure the acoustic phase 0(0 ~ 
for adiabatic models) and how it changes with scale just as the temperature but like the EE power spectrum [g^ has 
the added benefit of probing slightly larger scales than the first temperature peak. This property can help distinguish 
adiabatic and isocurvature models due to causal constraints on the generation of acoustic waves at the horizon at last 
scattering p8| . 

Finally, polarization also increases the viscosity of the fluid by a factor of 6/5, which has significant effects for the 
temperature. Even though the viscous imperfections of the fluid are small in the tight coupling region they can lead 
to significant dissipation of the fluctuations over time (see §[VD). 



C. Entropy and Heat Conduction 
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also represent imperfections in the fluid 



Differences in the bulk velocities of the photons and baryons Q( 
that lead to entropy generation and heat conduction. The baryon Euler equations (|68|) and (|7d) give 
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which may be iterated to the desired order in 1/f. For scalar fluctuations, this slippage leads to the generation of 
non-adiabatic pressure or entropy fluctuations 
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T^B = [Sp-yB - ct^B^PlB)IV-iB 



(96) 



as the local number density of baryons to photons changes. Equivalently, this can be viewed as heat conduction in 
the fluid. For vorticity fluctuations, these processes do not occur since there are no density, pressure, or temperature 
differentials in the fluid. 



D. Dissipation 



The generation of viscosity and heat conduction in the fluid dissipates fluctuations through the Euler equations 
with (|o|) and (||), 
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where we have dropped the aja factors under the assumption that the expansion can be neglected during the dissipation 
period. We have also employed Eqn. ( p6| ) to eliminate higher order terms in the vector equation. With the continuity 
equation for the scalars Gq — — fcO-j^ /3 — $ (see Eqn. ^ ^ = 0, to = ), we obtain 
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which is a damped forced oscillator equation. 

An interesting case to consider is the behavior in the absence of metric fluctuations ^, $, and V. The result, 
immediately apparent from Eqn. ( [9 7] ) and (|9§|), is that the acoustic amplitude {m = 0) and vorticity (to = 1) damp 
as exp[— (fc/fc]^')2] where 
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Notice that dissipation is less rapid for the vectors compared with the scalars once the fluid becomes baryon dominated 
i? ^ 1 because of the absence of heat conduction damping. In principle, this allows vectors to contribute more CMB 
anisotropics at small scales due to fluid contributions. In practice, the dissipative cut off scales are not very far apart 
since i? ^ 1 at recombination. 

Vectors may also dominate if there is a continual metric source. There is a competition between the metric source 
and dissipational sinks in Eqns. (p7|) and (Bq). Scalars retain contributions to 0q + ^ of ©(i?^, (^ — ^)/k'^) (see 



Eqn. (85) and p9[). The vector solution becomes 



Jo 
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which says that if variations in the metric are rapid compared with the damping then Q\ 
not occur. 



(100) 
V and damping does 



V. SCALING STRESS SEEDS 

Stress seeds provide an interesting example of the processes considered above by which scalar, vector, and tensor 
metric perturbations are generated and affect the temperature and polarization of the CMB. They are also the means 
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FIG. 7. Metric fluctuations from scaling anisotropic stress seeds sources. The same anisotropic stress seed (bold solid lines 
a^TVs (X fs/x) produce qualitatively different scalar (short-dashed), vector (long-dashed), and tensor (solid) metric perturba- 
tions. As discussed in the text the behavior scales with the characteristic time of the source Xc oc B^^ . The left panel (a) 
shows a source which begins to decay as soon as causally permitted (Bi = 1) and the right panel (b) the effect of delaying the 
decay (Bi = 0.2). We have displayed the results here for a photon-dominated universe for simplicity. 



by which cosmological defect models form structure in the universe. As part of the class of isocurvature models, all 
metric fluctuations, including the (scalar) curvature perturbation are absent in the initial conditions. To explore the 
basic properties of these processes, we employ simple examples of stress seeds under the restrictions they are causal 
and scale with the horizon length. Realistic defect models may be constructed by superimposing such simple sources 
in principle. 



We begin by discussing the form of the stress seeds themselves (§V A) and then trace the processes by which they 
form metric perturbations (§|VB|) and hence CMB anisotropics (§VC). 



A. Causal Anisotropic Stress 

Stress perturbations are fundamental to seeded models of structure formation because causality combined with 
energy-momentum conservation forbids perturbations in the energy or momentum density until matter has had the 
opportunity to move around inside the horizon (see e^. |30[| ). Isotropic stress, or pressure, only arises for scalar 
perturbations and have been considered in detail by |27[ . Anisotropic stress perturbations can also come in vector 
and tensor types and it is their effect that we wish to study here. Combined they cover the full range of possibilities 
available to causally seeded models such as defects. 

We impose two constraints on the anisotropic stress seeds: causality and scaling. Causality implies that correlations 
in the stresses must vanish outside the horizon. Anisotropic stresses represent spatial derivatives of the momentum 
density and hence vanish as fc^ for kr] <^1. Scaling requires that the fundamental scale is set by the current horizon 
so that evolutionary effects are a function oi x ~ krj. A convenient form that satisfies these criteria is ||27|,|3l| 



with 



fsi^) 



Ga^TT^") = 


--A(^)r"'iB{x), 


6 


sin(_Bix) sin(_B2a;) 
{Bix) {B2x) 


Bl - Bl 



(101) 



(102) 



with < (_Bi, B2) < 1. We caution the reader that though convenient and complete, this choice of basis is not optimal 
for representing the currently popular set of defect models. It suffices for our purposes here since we only wish to 
illustrate general properties of the anisotropy formation process. 
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FIG. 8. Temperature and polarization power spectra for a scaling anisotropic stress seeds with the minimal characteristic 
time Bi — 1 for scalars (S, solid), vectors (V, short-dashed), and tensors (T, long-dashed). Scalar temperature fluctuations at 
intermediate scales are dominated by acoustic contributions which then damp at small scales. B-parity polarization contribu- 
tions are absent for the scalars, larger by an order of magnitude than _E-parity contributions for the vectors and similar to but 
smaller than the _E-parity for the tensors. Features in the vector and tensor spectra are artifacts of our choice of source and 
are unlikely to be present in a realistic model. The background cosmology is set to ilo = 1, /i = 0.5, Qsh^ = 0.0125. 

Assuming Bi > B2, Bi controls the characteristic time after horizon crossing that the stresses are generated, i.e. 
the peak in Jb scales as krjc = Xc oc B^^ (see Fig. R). B2 controls the rate of decline of the source at late times. In 
the general case, the seed may be a smii of different pairs of {Bi,B2) which may also differ between scalar, vector, 
and tensor components. 

B. Metric Fluctuations 

Let us consi der h ow the anisotropic stress seed sources generate scalar, vector and tensor metric fluctuations. The 
form of Eqn. (102) implies that the metric perturbations also scale so that A:"^|/ip = f{x = krj) where / may be 
different functions for h = (^, $, V, H). Thus scaling in the defect field also implies scaling for the metric evolution 
and consequently the purely gravitational effects in the CMB as we shall see in the next section. Scattering introduces 
another fundamental scale, the horizon at last scattering r|^,, which we shall see breaks the scaling in the CMB. 

It is interesting to consider differences in the evolutions for the same anisotropic stress seed, A^™-' — 1 with Bi and 
B2 set equal for the scalars, vectors and tensors. The basic tendencies can be seen by considering the behavior at 



early times x 



< 



If X ^ 1 as well, then the contributions to the metric fluctuations scale as 



e'^V/JB = 0{x^) 



p/^^/fB - Oix-') 
e'^H/fB^O{x'), 



(103) 



where Jb — x^ for a; ^ 1. Note that the sources of the scalar fluctuations in this limit are the anisotropic stress 
and momentum density rather than energy density (see Eqn. pq ). This behavior is displayed in Fig. Wa). For the 
scalar and tensor evolution, the horizon scale enters in separately from the characteristic time Xc- For the scalars. 
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FIG. 9. Temperature-polarization cross power spectrum for the model of Fig. pi Independent of the nature of the sources, 
the cross power at angles larger than that subtended by the horizon at last scattering is negative for the scalars and vectors 
and positive for the tensors. The more complex structure for the scalars at small angular scales reflects the correlation between 
the acoustic efl^ective temperature and velocity at last scattering. 



the stresses move matter around and generate density fluctuations as ps 



The result is that the evolution of 



^ and $ steepens by x^ between 1 "^ x "^ Xc- For the tensors, the equation of motion takes the form of a damped 
driven oscillator and whose amplitude follows the source. Thus the tensor scaling becomes shallower in this regime. 
For x^ Xc both the source and the metric fluctuations decay. Thus the maximum metric fluctuation scales as 



k^'^V/fB{Xc) 



0{x\)., 
0{xl)., 



k"^^/!B{xc)^0{x\), 
k'/^HlfBixc) = 0{x^^). 



(104) 



For a late characteristic time a;c > 1, fluctuations in the scalars are larger than vectors or tensors for the same source 
(see Fig. CTj). The ratio of acoustic to gravitational redshift contributions from the scalars scale as x~^ by virtue of 
pressure support in Eqn. ([Sq) and thus acoustic oscillations become subdominant as Bi decreases. 



C. CMB Anisotropies 



Anisotropy and structure formation in causally seeded, or in fact any isocurvature model, proceeds by a qualitatively 
different route than the conventional adiabatic inflationary picture. As we have seen, fluctuations in the metric are 
only generated inside the horizon rather than at the initial conditions (see §V_B). Since CMB anisotropies probe 
scales outside the horizon at last scattering, one would hope that this striking difference can be seen in the CMB. 
Unfortunately, gravitational redshifts between last scattering and today masks the signature in the temperature 
anisotropy. The scaling ansatz for the sources described in § V A in fact leads to near scale invariance in the large angle 
temperature because fluctuations are stimulated in the same way for each A:-mode as it crosses the horizon between 
last scattering and the present. While these models generically leave a different signature in modes which cross the 
horizon before last scattering p3,E^, models which mimic adiabatic inflationary predictions can be constructed |31|. 

Polarization provides a more direct test in that it can only be generated through scattering. The large angle 
polarization reflects fluctuations near the horizon at last scattering and so may provide a direct window on such 
causal, non-inflationary models of structure formation. One must be careful however to separate scalar, vector and 
tensor modes whose different large angle behaviors may obscure the issue. Let us now illustrate these considerations 
with the specific examples introduced in the last section. 

The metric fluctuations produced by the seed sources generate CMB anisotropies through the Boltzmann equation 
(|60|). We display an example with Bi = 1 and B2 = 0.5 in Fig. |8[ Notice that scahng in the sources does indeed lead 
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FIG. 10. Same as Fig. N except with a larger characteristic time Bi — 0.2, B2 — 0.1. Scalar gravitational redshift effects now 
dominate over scalar acoustic as well as vector and tensor contributions for the same stress source due the process by which 
stress perturbations generate metric fluctuations (see Fig. |7|). 



to near scale invariance in the large angle temperature but not the large angle polarization. The small rise toward the 
quadrupole for the tensor temperature is due to the contribution of long-wavelength gravity waves that are currently 
being generated and depends on how rapidly they are generated after horizon crossing. Inside the horizon at last 
scattering (here £ > 200), scalar fluctuations generate acoustic waves as discussed in §|rVj which dominate for small 
characteristic times a;c ~ 1. On the other hand, these contributions are strongly damped below the thickness of the 
last scattering surface by dissipational processes. Note that features in the vector and tensor spectrum shown here are 
artifacts of our choice of source function. In a realistic model, the superposition of many sources of this type will wash 
out such features. The general tendencies however do not depend on the detailed form of the source. Note that vector 
and tensor contributions damp more slowly and hence may contribute significantly to the small-angle temperature 
anisotropy. 

Polarization can only be generated by scattering of a quadrupole temperature anisotropy. For seeded models, 
scales outside the horizon at last scattering fcry* ^ 1 have not formed significant metric fluctuations (see e.g. Fig. m. 
Hence quadrupole fluctuations, generated from the metric fluctuations through Eqns. (^l|), (|9^), and (p3), are also 
suppressed. The power in k of the polarization thus drops sharply below fcry* = 1. This drop of course corresponds to a 
lack of large angle power in the polarization. However its form at low £ depends on geometric aspects of the projection 
from k to £. In these models, the large angle polarization is dominated by projection aliasing of power from small 
scales fcry* ^ 1. The asymptotic expressions of Eqn. (pq) thus determine the large angle behavior of the polarization: 
I'^Ce oc £^ for scalars (EE), £'^ for vectors and £'^ for tensors [EE and BB)\ the cross spectrum [QE) goes as £'^ for 
each contribution. For comparison, the scale invariant adiabatic inflationary prediction has scalar polarization {EE) 
dropping off as £'^ and cross spectrum (9-E) as £'^ from Eqn. (|9^) because of the constant potential above the horizon 
(see Fig. |6|) . Seeded models thus predict a more rapid reduction in the scalar polarization for the same background 
cosmology. 

Polarization can also help separate the three types of fluctuations. In accord with the general prediction (see § |III D| ) , 
scalars produce no i?-parity polarization, whereas vector _B-parity dominates -E-parity polarization by a factor of 6 
and tensor i?-parity is suppressed by a factor of 8/13 (see Eqn. d23)). Differences also arise in the temperature- 
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polarization cross power spectra Cf^ shown in Fig. g. Independent of the nature of the source, above the angle 
the horizon subtends at last scattering, scalar and vector temperature perturbations from the last scattering surface 



[ pi[ are anticorrelated with polarization, whereas they are correlated for tensor perturbations (see § |[VB| and [g|). 
Inside the horizon, the scalar polarization follows the scalar velocity which is tt/2 out of phase with the effective 
temperature (see Eqn. ^^. In the adiabatic model, scalar cross correlation reverses signs before the first acoustic 
peak, as compression overcomes the gravitational redshift of the Sachs- Wolfe effect, unlike the isocurvature models 
(see Fig. and pq]). The sign test to distinguish scalars from tensors must thus be performed on scales larger than 
twice the first peak. Conversely, to use the cross correlation to distinguish adiabatic from isocurvature fluctuations, 
the scalar and tensor contributions must be separated. 

How do these results change with the model for the seeds? As we increase the characteristic time Xc by decreasing 



Bi = 0.2, the main effect comes from differences in the generation of metric fluctuations discussed in §VB. For the 



same amplitude anisotropic stress, scalars contributions dominate the vector and tensor contributions by factors of 



B^ (see Eqn. 104 ) . Note however that the scalar contributions come from the gravitational redshifts between 



last scattering and today rather than the acoustic oscillations (see § [V A ) and hence produce no strong features. 
Because of the late generation of metric fluctuations in these models, the peak in the polarization spectra is also 
shifted with Xc- Note however that the qualitative behavior of the polarization described above remains the same. 

Although these examples do not exhaust the full range of possibilities for scaling seeded models, the general 
behavior is representative. Equal amplitude anisotropic stress sources tend to produce similar large angle temperature 
anisotropics if the source is active as soon as causally allowed Xc ~ I. Large angle scalar polarization is reduced as 
compared with adiabatic inflationary models because of causal constraints on their formation. This behavior is not 
as marked in vectors and tensors due to the projection geometry but the relative amplitudes of the i?-parity and 
_B-parity polarization as well as the QE cross correlation can be used to separate them independently of assumptions 
for the seed sources. Of course in practice these tests at large angles will be difficult to apply due to the smallness of 
the expected signal. 

Reionization increases the large angle polarization signal because the quadrupole anisotropics that generate it can 



be much larger |32|. This occurs since decoupling occurs gradually and the scattering is no longer rapid enough to 
suppress anisotropics. The prospects for separating the scalars, vectors and tensors based on polarization consequently 
also improve [ p3[ . 

For angles smaller than that subtended by the horizon at last scattering, the relative contributions of these effects 
depends on a competition between scalar gravitational and acoustic effects and the differences in the generation and 
damping behavior of the scalar, vector and tensor perturbations. 

VI. DISCUSSION 

We have provided a new technique for the study of temperature and polarization anisotropy formation in the CMB 
which introduces a simple and systematic representation for their angular distributions. The main virtue of this 
approach is that the gravitational and scattering sources are directly related to observable properties in the CMB. 
One can then explore properties that are independent of the source, which tell us the broad framework, e.g. the 
classical cosmological parameters and the nature of fluctuations in the early universe, and identify properties that 
are dependent on the source, which can help pin down the model for structure formation. An example of the former 
is the fact that scalar fluctuations generate no magnetic parity polarization |^,y, vectors generate mainly magnetic 
parity polarization, and tensors generate comparable but somewhat smaller magnetic parity polarization. Large angle 
polarization of the three components are also constrained by model-independent geometric arguments in its slope and 
its correlation with the temperature anisotropy. If the scalar contributions can be isolated from the vectors, tensors 
and other foreground sources of polarization from these and other means, these constraints translate into a robust 
distinction between isocurvature and adiabatic models for structure formation. 

In our representation, the temperature and polarization distributions are projections on the sky of four simple 
sources, the metric fluctuation (via the gravitational redshift), the intrinsic temperature at last scattering, the baryon 
velocity at last scattering (via the Doppler effect) and the temperature and polarization quadrupoles at last scattering 
(via the angular dependence of Compton scattering). As such, it better reveals the power of the CMB to probe 
the nature of these sources and extract information on the process of structure formation in the universe. As an 
example, we have explored how general properties of scaling stress seeds found in cosmological defect models manifest 
themselves in the temperature and polarization power spectra. The framework we have provided here should be useful 
for determining the robust signatures of specific models for structure formation as well as the reconstruction of the 
true model for structure formation from the data as it becomes available. 
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TAB-2. Commonly used symbols, m = 0, ±1, ±2 for the scalars, vectors, and tensors. For the fluid variables / — > 7 
for the photons, / — > i? for the baryons and ^B for the photon-baryon fluid. X = Q, E, B for the temperature- 
polarization power spectra. 
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